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Abstract 

We review the use of Legendre transforms in the formulation of electrostatic energies in 
condensed matter. We show how to render standard functionals expressed in terms of the 
electrostatic potential, cj), convex - at the cost of expressing them in terms of the vector 
field D. This leads to great simplification in the formulation of numerical minimisation of 
electrostatic energies coupled to other physical degrees of freedom. We also demonstrate 
the equivalence of recent functionals for dielectrics derived using field theory methods to 
classical formulations in terms of the electric polarisation. 

I. Introduction 

The Legendre transform is a powerful tool with multiple applications in physics [T] . In classical 
mechanics it allows one to interchange the Lagrangian and Hamiltonian viewpoints; in ther- 
modynamics one regularly transforms ensembles to simplify calculations, choosing the ensemble 
which most closely idealizes a given experimental setup. In this article we will demonstrate the 
utility of Legendre transforms in reformulating energies and free energies in electrostatics, with 
a particular eye for numerical applications and mean field theory. 

Our principal motivation for a deeper study of this transformation applied to electrostatic 
problems is quite practical. Many formulations of (free) energy functions in condensed matter 
physics involve the electrostatic potential, an important example is the Poisson-Boltzmann energy 
functional in the theory of ionic solutions. However when we examine closely these functionals 
we see that they are concave functions of the potential. While the stationary value of the 
functional is indeed the correct value of the energy that we wish to study, the concavity leads to 
complications in many situations. In particular we can not perform a simultaneous minimisation 
of both electrostatic and configurational energies in a simulation. One is generally obliged to 
fully solve by iterative methods the electrostatic problem at each time step of the iteration over 
the configurational degrees of freedom - such as densities or polymer configuration. This leads 
to codes which are complicated to write, and sometimes slow to run. 

We remind the reader that for a convex function f(x), its Legendre transform is defined [T] 
from the expression 

C[f](s)=g(s) = sx-f(x) (1) 

where on the right of eq. ([!} we express x as a function of s from the equation s = This 
transformation is an involution: C [g] (x) = f(x). The simplest example is an Hookian spring for 
which f(x) = kx 2 /2 ; then the transformation of eq. ([T]) gives g(s) = s 2 /(2k). We will also use 
the notation £(/) = /. In this article we show that by introducing new variational parameters 
in a free energy with the help of Lagrange multipliers and then performing a Legendre transform 
of the resulting free energy we can find functionals that are convex in all degrees of freedom; 
we will illustrate this with a mean field formulation of phase separation coupled to electrostatic 
interactions. 

We now illustrate the approach with the simplest possible electrostatic problem, interaction 
between free charges, pf in a heterogeneous dielectric medium: Consider the energy functional 
expressed in terms of the electric potential <p. 
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The variational equation for the field is then the Poisson equation: 

div e(r)grad 4> = —Pfi r ) (3) 
Substituting the solution of the Poisson equation in the electrostatic energy we find 

U=\J^(r)p f (r)d 3 r (4) 

We convert the variational problem for the potential by introducing the new variable E = — V0. 
To do so we introduce the (vector) Lagrange multiplier D. The stationary point of eq. ([2J is 
identical to that of the following expression [2J: 

C/ = ||-^^+P/0 + D.(E + V0)| d 3 r (5) 

It is at this point that we recognise that the variational equations for E correspond to a Legendre 
transform with dual variable D. We also integrate by parts the product D • V(j> to find — 0div D, 
dropping boundary terms assumed to be zero. Thus the stationary point of eq. ([2]) is identical 
to the stationary point of 

U = J j-^y + 0(p-divD)j d 3 r (6) 

Variations in <p now impose Gauss' law, div D — pf = 0, while the energy has been rendered 
convex by the transformations introduced, [3J 0] ■ 

We now illustrate applications of this transformation to two problems: Firstly the theory 
of phase separation of immiscible fluids in the presence of electrostatic interactions due to ions. 
Secondly we provide a translation between two very different visions of the theory of dielectrics. 
Recent formulations of implicit dielectrics pass by elaborate field-theory mappings and find a 
generalised Poisson-Boltzmann equation with a Langevin correction. We show how to map this 
description onto a free energy expressed in terms of a polarisation field with long-ranged dipolar 
interactions. We believe that these equivalent descriptions can lead to a deeper understanding 
of the underlying physics. 

In the following we will work with free energy densities, rather than the integrated energies 
and we (silently) integrate by parts when needed. 



II. Phase separation coupled to electrostatics 

A mixture of two solvents (A and B) near their miscibility limit and in the presence of salt displays 
interesting properties which have been explored in recent experiments [S]- Density fluctuations 
couple to the dielectric properties of the medium, and in turn influence the partition of ions in 
the fluctuating solvent field. The experimental system has turned out to be very rich, and allows 
one to adjust the effective interaction between colloidal particles using temperature as a control 
parameter. 

A simplified theoretical description of such systems is given in [5J[7] who propose the following 
free energy density expressed in terms of the densities and potentials: 

/(^,c+,c_) = - ie(*)(V^) 2 + (c+ -c_)e0 

- (Au+c + + Au- C _) * + k B T^2( Cj ln( Cj /c j0 ) - cj) (7) 

3 

describes the composition fluctuations of the fluid mixture. c+ and c_ arc the concentration 
of positive and negative monovalent ions, with Au + and Au~ their relative preferences between 
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a A-liquid environment and a B-liquid environment. As above <f> is the electrostatic potential. 
We see that fluctuations in concentration couple via e( v E') to a coupling with the concentration 
fluctuations of the ions. 

/ m (*) includes all the terms that are only dependent on /,„(*) = /o(*) + f (V*) 2 - 
with fo(^) the free energy due to the mixing of the two solvents. It can, for exemple, be written 
as a binary mixture free energy density: /o(^) oc 'J/log^) + (1 — 4 r )log(l — + x*(l — 
with x the Flory parameter [HIE], or as a Landau expansion fo(^) oc a — ^c) 2 + 7 — ^c) 4 , 
with a being temperature dependent, 7 positive, and $ c the critical composition [6]. 

Optimising eq. (0 over c+ and c_ , the density becomes : 

/(*, 0) = / m (*)-ie(*)(V0) 2 - fc s Tc 0+ exp(/3A M +* - pe<f>) 

-k B Tco-exp(PAu-^f + Pecfy (8) 

With a symmetric electrolyte: co+ = cq— , and if we assume the ions are similar in their 
interaction with the solvents : Au + = Ait~, /(\I/,<^) simplifies into : 

0) = / w (¥) - I e (*)(V0) 2 - 2k B Tc Q exp(pAu^) cosh(/3e0) (9) 

We recognise here a generalisation of the well known Poisson-Boltzmann functional for a symmet- 
ric electrolyte. The description is adapted to analytical solutions but in the monophase region 
of the phase diagram <f>) is convex in \Jj but concave in <f>. In complicated geometries if one 
wishes to minimize this free energy numerically one has to solve saddle point equations, simple 
minimization will not give the correct answer. We now implement the transformation introduced 
above from the potential <j> to the electric displacement D and use the fact that the Legendre 
transform of cosh is : 



C[Acosh(B<t>)](Q = A 

=Ag H/AB) 

After some calculation we find 



f/(AB)asinh(£/AB) - \J(Z/AB) 2 + 1 



(10) 



/(*, D) = /„(*) + + 2k B Tc Q e^ 9 ^^-^ j (11) 

We have thus reached our objective: we have built an equivalent description of the system with 
the stationary conditions conserved and a local and convex function. The disadvantage is that 
there are more degrees of freedom in the vector field D than in the scalar field <f>, but the 
advantage is that a global minimising principle can be used and the functional can be directly 
programmed for the solution of the coupled electrostatic-phase separation problem. 

We note that mean field description of the packing of DNA in a virus [TU] contains many 
similar theoretical features and is also amenable to similar transformations. In this problem the 
field W corresponds to the square root of the monomer density. 

III. From Poisson-Langevin to Polarization 

We now consider theories of explicit Langevin dipoles and how these can be incorporated into 
the formulation of the free energy in terms of convex free energy functions. Recent work on 
improving the description of solvation of proteins [11] has considered an explicit model for the 
solvent in terms of Langevin dipoles. If we neglect the volume of ions and dipoles they find the 
free energy density for a mixture of symmetric ions and neutral dipoles : 

6 O (V0) 2 sinh(/? Po |V0|) 

/ = Pf<p 2 lon (P q w " dip — pp \v^)\ — ' 
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where to simplify the presentation we have neglected effects of finite ion size. The parameters A 
are related to the chemical activities of the ions and the dipoles. 

As in previous work |12j we start by using a Lagrangian multiplier, D to replace (V</>) by its 
electrostatic equivalent — E. We find 

/ = P/*- f y--ff(0-fc(E) + D.(V^ + E) (13) 

where h(E) is the free energy density due to the dipoles and g((j>) the free energy due to free 
ions. We now diverge from our previous treatment and introduce a new variable P which we 
will show is the physical polarisation variable. We do this by performing a Legendre transform 
on h(E) to find h(P); we then find 

/ = 4>{pj - div D) - - gtf) + h(P) + E • (D - P) (14) 

Clearly by definition of the Legendre transform performing variations with respect to P on 
eq. mijl gives eq. (1T31) . 

We now perform two more transforms firstly to eliminate the potential, but secondly to 
eliminate the electric field E. We find 

(D-P) 2 ~ 

/ = 1 2 ' + MP) + 9(Pf - div D) (15) 

In the absence of free ions the function g reduces to the constraint of Gauss' law. This is exactly 
the form postulated in |13j . It is particularly transparent for understanding the physical limits 
on response functions [141 115) and the origin of the negative dielectric constant observed in 
structured fluids. 

We will now work from eq. (115|) to demonstrate its equivalence to other formulations of 
electrostatic interactions expressed in terms of the polarization P. To do this we will eliminate 
the variable D, which will bring us back to other more familiar forms for the electrostatic energy 
at the cost of re-introducing long-ranged dipolc-dipole interactions between the polarization 
variables. 



Eliminating the Displacement field 

Let us work in the limit where linear response is valid, in which case we expand h to quadratic 
order: 

MP) = ^— (16) 

where \ 1S a material parameter. Taking variations of eq. (|14j) with respect to P and then E we 
find that 

P = e xE, e E = D P (17) 

Thus the parameter x is t ne electric susceptibility of the medium. The polarisation variable 
is indeed playing the role we expect from standard treatments of Maxwell's equations. The 
free energy of the fluctuating dipoles (in the absence of free ions) can then be found from the 
functional 

(D-P) 2 P 2 
^S^ + 2^)-^ divD -^ < 18 > 

where the last term is a Lagrange multiplier for the constraint of Gauss' law. On taking variations 
of eq. (fT8f with respect to D we find that 

D-P = -e o V0 (19) 
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Thus the free energy density can be written as 

(Vd>) 2 P 2 

f = eo { -^- + ^— (20) 

where 

e V 2 (j) = -p f + div P 

0(r) = / 47r eo |r-r'| W) ^ ^ ^ 

We now substitute eq. (f2~Tj) in eq. (|20l) and introduce the bare electric field Eo as follows: Eo = 
— V0o with eoV 2 0o = — Pf- The free energy density can then be expressed in terms of the 
polarization as 

1 fdivP(r)divP(r') - - , 



U = - — ^ d A Y d?r 

2 J 4^e |r-r'| 

/(!2M_E -P + ^-l d 3 r (22) 
7 12 2e oX (r) / V ; 

This formulation of the free energy is widely used in theoretical chemistry and is that used by 

m 

The energy eq. (|22p can be expressed in an even more physically transparent manner by 
integrating by parts the first double integral to transfer the derivatives from the polarization to 
the function l/\r — r'|. If we do this we find that the double integral is transformed to 

i f P(r)T(r - r')P(r') d 3 r d 3 r (23) 

where the dipole operator is 



T(r) 



l-3|r)(r| 5(r) 



(24) 



Anr 3 eo 3eo 

We now proceed in a more abstract manner, considering that the polarisation variables are 
assembled into a vector and the dipolar interactions form a matrix, T. Eq. (1221) is simply a 
quadratic form in P, so that 

P(K + T)P En , x 

U= [ - 2 ^ -E .P + £ o^ (25) 

with the diagonal matrix eoK r r = Xr 1 - If we now calculate the response of the polarization field 
to an external perturbation Eo we find 

P = - 1 - E (26) 
K + T ° V ' 

The total electric field is then given by two contributions, the original imposed field Eo and 
that due to the dipole density P: 

E = E - TP (27) 
K 



K + T 

We also use 



E (28) 



D — P = e E (29) 
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to find 

D = ^£ Eo (30) 
K+T ' V ' 

Again all these equations are non-local - since they involve the long-ranged operator T but we 
find that 

D = (l + Xr )e E = eE (31) 

a purely local constitutive equation between the electric field and the electric displacement. 

We conclude that by careful study of the mean field equations coming from the formulation 
of the dielectric properties of a medium in eq. (|12p we have been able to derive the equivalence 
to the standard continuum formulation of electrostatic arising from Maxwell's equations. 



IV. Conclusions 

We have shown that the Legendre transform can be used to translate between multiple forms of 
the energy in mean-field theories. All the formulations are numerically equivalent but different 
forms put the emphasis on different degrees of freedom in electromagnetism. For numerical work 
it is advantageous to work with a formulation which is both convex and local. This is achieved in 
Poisson-Boltzmann theory by choosing the electric displacement D as the fundamental thermo- 
dynamic field. In this way all physical degrees of freedom can be treated in an equivalent manner 
in numerical solvers. It is no longer necessary to completely solve the electrostatic problem for 
each iteration of other external degrees of freedom. Very similar conclusions have also been found 
in quantum chemistry [T7] 

We have also demonstrated that energy functionals for dielectrics can be translated into 
equivalent forms by introducing the physical polarization. We then mapped a linearised form of 
the theory to the Marcus energy function, widely used in the theoretical chemistry literature. 
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